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Rate-dependent plasticity: creep and swelling 


The classical deviatoric metal creep behavior in Abaqus/Standard: 


can be defined using user subroutine CREEP or by providing parameters as input 
for some simple creep laws; 


can model either isotropic creep (using Mises stress potential) or anisotropic creep 
(using Hill's anisotropic stress potential); 


is active only during steps using the coupled temperature-displacement procedure, 
the transient soils consolidation procedure, and the quasi-static procedure; 


requires that the material's elasticity be defined as linear elastic behavior; 


can be modified to implement the auxiliary creep hardening rules specified in 
Nuclear Standard NEF 9-5T, “Guidelines and Procedures for Design of Class 1 
Elevated Temperature Nuclear System Components”; these rules are exercised by 
means of a constitutive model developed by Oak Ridge National Laboratory (ORNL 
- Oak Ridge National Laboratory constitutive model); 


can be used in combination with creep strain rate control in analyses in which the 
creep strain rate must be kept within a certain range; and 

can potentially result in errors in calculated creep strains if anisotropic creep and 
plasticity occur simultaneously (discussed below). 


Rate-dependent gasket behavior in Abaqus/Standard: 


uses unidirectional creep as part of the model of the gasket's thickness-direction 
behavior; 


can be defined using user subroutine CREEP or by providing parameters as input 
for some simple creep laws; 


is active only during steps using the quasi-static procedure; and 


requires that an elastic-plastic model be used to define the rate-independent part 
of the thickness-direction behavior of the gasket. 


Volumetric swelling behavior in Abaqus/Standard: 


can be defined using user subroutine CREEP or by providing tabular input; 
can be either isotropic or anisotropic; 


is active only during steps using the coupled temperature-displacement procedure, 
the transient soils consolidation procedure, and the quasi-static procedure; and 


requires that the material's elasticity be defined as linear elastic behavior. 


The following topics are discussed: 


Creep behavior 

Volumetric swelling behavior 

User subroutine CREEP 

Removing creep effects in an analysis step 
Integration 

Loading control using creep strain rate 
Elements 

Output 


Products: Abaqus/Standard Abaqus/CAE 


Creep behavior 


Creep behavior is specified by the equivalent uniaxial behavior—the creep “law.” In practical cases creep laws are typically of 
very complex form to fit experimental data; therefore, the laws are defined with user subroutine CREEP, as discussed below. 
Alternatively, five common creep laws are provided in Abaqus/Standard: the power law, the hyperbolic-sine law, the double 
power law, the Anand law, and the Darveaux law. These standard creep laws are used for modeling secondary or steady-state 
creep. Creep is defined by including creep behavior in the material model definition (Material data definition). Alternatively, 
creep can be defined in conjunction with gasket behavior to define the rate-dependent behavior of a gasket. 


Input File Usage: 
Use the following options to include creep behavior in the material model definition: 


*MATERIAL 
* CREEP 


Use the following options to define creep in conjunction with gasket behavior: 
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See Also 
About the material library 
Inelastic behavior 


In Other Guides 


Defining the gasket behavior 
directly using a gasket behavior 
model 


*CREEP 


*CREEP STRAIN RATE 
CONTROL 


*POTENTIAL 
*SWELLING 
*RATIOS 

Defining a creep law 
Defining swelling 
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*GASKET BEHAVIOR 
* CREEP 


Abaqus/CAE Usage: 
Property module: material editor: Mechanical > Plasticity > Creep 
Choosing a creep model 


The power-law creep model is attractive for its simplicity. However, it is limited in its range of application. The time-hardening 
version of the power-law creep model is typically recommended only in cases when the stress state remains essentially 
constant. The strain-hardening version of power-law creep should be used when the stress state varies during an analysis. In 
the case where the stress is constant and there are no temperature and/or field dependencies, the time-hardening and strain- 
hardening versions of the power-creep law are equivalent. For either version of the power law, the stresses should be 
relatively low. 


In regions of high stress, such as around a crack tip, the creep strain rates frequently show an exponential dependence of 
stress. The hyperbolic-sine creep law shows exponential dependence on the stress, 0, at high stress levels (c/o° > 1, where 


o° is the yield stress) and reduces to the power-law at low stress levels (with no explicit time dependence). 


The double power, Anand, and Darveaux models are particularly well suited for modeling the behavior of solder alloys used in 
electronic packaging and have been shown to produce accurate results for a wide range of temperatures and strain rates. 


None of the above models is suitable for modeling creep under cyclic loading. The ORNL model (ORNL - Oak Ridge National 
Laboratory constitutive model) is an empirical model for stainless steel that gives approximate results for cyclic loading 
without having to perform the cyclic loading numerically. Generally, creep models for cyclic loading are complicated and must 
be added to a model with user subroutine CREEP or with user subroutine UMAT. 


Modeling simultaneous creep and plasticity 


If creep and plasticity occur simultaneously and implicit creep integration is in effect, both behaviors may interact and a 
coupled system of constitutive equations needs to be solved. If creep and plasticity are isotropic, Abaqus/Standard properly 
takes into account such coupled behavior, even if the elasticity is anisotropic. However, if creep and plasticity are anisotropic, 
Abaqus/Standard integrates the creep equations without taking plasticity into account, which may lead to substantial errors in 
the creep strains. This situation develops only if plasticity and creep are active at the same time, such as would occur during a 
long-term load increase; one would not expect to have a problem if there is a short-term preloading phase in which plasticity 
dominates, followed by a creeping phase in which no further yielding occurs. Integration of the creep laws and rate-dependent 
plasticity are discussed in Rate-dependent metal plasticity (creep). 


Time power law and power law models 


The time power law and power law models described below are equivalent to the “time hardening” and the “strain hardening” 
forms but avoid their drawbacks. The time power law and power law models rewrite the laws in such a way that the typical 
parameter values do not cause numerical difficulties. In addition, the units of all of the parameters are physical, which makes 
unit conversion easier if it is required. 


Time power law model 


The time power law model has the following form: 


= RH in 

+Cr i q . 

= Eq] — Egt |] , 
qo 


Cr ~ . 
where € and q are defined above; and €g, n, m, and qo are material parameters. 


The model is equivalent to the “time hardening” form. It is recommended that you use the time power law model when the 


value of the parameter A is very small (A < 10°?" ), In this case the equivalent time power law model is obtained by setting 
m+1 1 


do = £g” Ae, keeping the parameters n and ™m unchanged, and setting €9 to an arbitrary value greater than zero 
(typically, €9 is set to one). 


Input File Usage: 





*CREEP, LAW=TI POWER LAW 


Abaqus/CAE Usage: 











Property module: material editor: Mechanical > Plasticity > Creep: Law: Time Power 
Power law model 
The power law model has the following form: 


1 
m+1 


. q 
—E =é9|/|— (m+ 1)E7 : 
qo 





~ n m 
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CP on . 
wheree , Ev, and q are defined above; and €9, n, m, and qo are material parameters. 


This model is equivalent to the “strain hardening” form. It is recommended that you use the power law model when the value 


of the parameter A is very small (A < 10-7"), In this case the equivalent power law model is obtained by setting 
picasa 1 
go = €9” A=, keeping the parameters n and m unchanged, and setting €9 to an arbitrary value greater than zero 


(typically, Eg is set to one). 
Input File Usage: 
*CREEP, LAW=POWER LAW 
Abaqus/CAE Usage: 
Property module: material editor: Mechanical > Plasticity > Creep: Law: Power 


Time/ Strain hardening models 


Time hardening and strain hardening models to specify creep are available. However, to avoid the drawbacks of these models, 
it is recommended that you use the time power law and power law models (see Time power law and power law models). 


Time hardening form 


The “time hardening” form is the simpler of the two forms and is defined as 


cr ~ 
—E = Aq”t™, 
where 
cr 
E 
j ee : F 2 -cr scr 
is the uniaxial equivalent creep strain rate, ze 7€, 
q 
is the uniaxial equivalent deviatoric stress, 
t 


is the total or the creep time, and 
A, n, and m 
are defined by you as functions of temperature. 
q is Mises equivalent stress or Hill's anisotropic equivalent deviatoric stress according to whether isotropic or anisotropic creep 
behavior is defined (discussed below). For physically reasonable behavior A and n must be positive and -—1 <<m< 0. 
Input File Usage: 
*CREEP, LAW=TIME 
Abaqus/CAE Usage: 
Property module: material editor: Mechanical > Plasticity > Creep: Law: Time-Hardening 
Strain hardening form 


The “strain hardening” form is 


1 


Ee" = (Ag"[(m + 1)E7]") ™, 


where €~ and gq are defined above and E“ is the equivalent creep strain. 
Input File Usage: 
*CREEP, LAW=STRAIN 
Abaqus/CAE Usage: 
Property module: material editor: Mechanical > Plasticity > Creep: Law: Strain-Hardening 
Numerical difficulties 


Depending on the choice of units for either form, the value of A may be very small for typical creep strain rates. If A is less 
than 10727, numerical difficulties can cause errors in the material calculations. Therefore, use another system of units or use 
the time power law or power law model (described below) to avoid such difficulties in the calculation of creep strain 
increments. 


Time-dependent behavior 


In the “time hardening” form and the time power law model, the total time or the creep time can be used. The total time is 
the accumulated time over all general analysis steps. The creep time is the sum of the times of the procedures with time- 
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dependent material behavior. If the total time is used, it is recommended that small step times compared to the creep time be 
used for any steps for which creep is not active in an analysis; this is necessary to avoid changes in hardening behavior in 
subsequent steps. 


Input File Usage: 
Use one of the following options: 


*CREEP, TIME=TOTAL (default) 
*CREEP, TIME=CREEP 


Abaqus/CAE Usage: 














Property module: material editor: Mechanical > Plasticity > Creep: Time: Total or Creep 
Hyperbolic-sine law model 


The hyperbolic-sine law is available in the form 


Cr n AH 
é = A(sinh Bq)” exp| -————— }], 
R(0—- 64) 


where 
é" and q 
are defined above, 
0 
is the temperature, 
02 
is the user-defined value of absolute zero on the temperature scale used, 
AH 
is the activation energy, 
R 
is the universal gas constant, and 
A, B, andn 
are other material parameters. 


This model includes temperature dependence, which is apparent in the above expression; however, the parameters A, B, n, 
A\H, and R cannot be defined as functions of temperature. 


Input File Usage: 
Use both of the following options: 


*CREEP, LAW=HYPERB 
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=67 


Abaqus/CAE Usage: 
Define both of the following: 
Property module: material editor: Mechanical > Plasticity > Creep: Law: Hyperbolic-Sine 
Any module: Model > Edit Attributes > model_name: Absolute zero temperature 


Numerical difficulties 


As with the power law, A may be very small for typical creep strain rates. If A is very small (such as less than 10-27), use 
another system of units to avoid numerical difficulties in the calculation of creep strain increments. 


Anand model 


The Anand model is available in the form 


* = alsinn(¢2)]" exp ww , 


where 
cr 


E ,q,R,9, and 67 
are defined above, 


Q 


is the activation energy, 
S 


https://help.3ds.com/2020/english/dssimulia_established/simacaematrefmap/simamat-c-ratedepcreep.htm?contextscope=all&format=P&value=undefi... 4/11 


1/3/22, 4:21 PM Rate-dependent plasticity: creep and swelling - 2020 - dssimulia_established 


is the deformation resistance, and 
A,m, and & 


are material parameters. 


The evolution equation for the deformation resistance, s (initially s = Sq), is 





with 


’] 


where 
ho = Ay + A; (0 — 67) + An (0 — 67)” + Agé™ + Ay(E)’, 


anda, n, 8, Ag, Ai, Az, A3, and Ag are material parameters. 


In addition, the initial deformation resistance is a function of temperature of the form 
2 
59 = $1 + Sp (0-07) + S3(0- 67) , 


where $1, S2, and $3 are material parameters. 


Input File Usage: 
Use both of the following options: 


*CREEP, LAW=ANAND 
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=07 


Abaqus/CAE Usage: 
Define both of the following: 
Property module: material editor: Mechanical > Plasticity > Creep: Law: Anand 
Any module: Model > Edit Attributes > model_name: Absolute zero temperature 
Darveaux model 


The Darveau model involves both primary and secondary creep. The secondary creep (steady-state) component is described 
by a standard hyperbolic sine law 


Q 
R(0— 07) 


cr 
Ss 


= C,,[sinh(aq)|” exp| — 


The steady-state law is modified to include the primary creep effects through 
ev = é, [1 + epBexp (—Bé; t)] ; 


where 

Ee, G, R, Q, 8, and 67 
are defined above, 

Css 


is the steady-state creep prefactor, 
a 


is the steady-state creep power law breakdown, and 
n,€er,and B 
are other material parameters. 
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Input File Usage: 
Use both of the following options: 


*CREEP, LAW=DARVEAUX 
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=07 


Abaqus/CAE Usage: 
Define both of the following: 
Property module: material editor: Mechanical > Plasticity > Creep: Law: Darveaux 
Any module: Model > Edit Attributes > model_name: Absolute zero temperature 
Double power model 


The double power law is available in the form 


ol B, q\"% By q\° 
“= Aeo(—G=ay) (zy) tAer(-gaay) (a) 


where 
Cr. 


— ,q,9, and 67 


are defined above, 
00 


is the normalized stress, and 
Al, A», Bi, Bo, Ci, and Cy 
are other material parameters. 
Input File Usage: 
Use both of the following options: 


*CREEP, LAW=DOUBLE POWER 
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=07 


Abaqus/CAE Usage: 
Define both of the following: 
Property module: material editor: Mechanical > Plasticity > Creep: Law: Double Power 
Any module: Model > Edit Attributes > model_name: Absolute zero temperature 
Anisotropic creep 


Anisotropic creep can be defined to specify the stress ratios that appear in Hill's function. You must define the ratios Ri; in 


each direction that will be used to scale the stress value when the creep strain rate is calculated. The ratios can be defined as 
constant or dependent on temperature and other predefined field variables. The ratios are defined with respect to the user- 
defined local material directions or the default directions (see Orientations). Further details are provided in Hill anisotropic 
yield/creep. Anisotropic creep is not available when creep is used to define a rate-dependent gasket behavior since only the 
gasket thickness-direction behavior can have rate-dependent behavior. 


Input File Usage: 


* POTENTIAL 
Abaqus/CAE Usage: 


Property module: material editor: Mechanical > Plasticity > Creep: Suboptions > Potential 


ry 
Volumetric swelling behavior 


As with the creep laws, volumetric swelling laws are usually complex and are most conveniently specified in user subroutine 
CREEP as discussed below. However, a means of tabular input is also provided for the form 


ee =F (Oi Fis farese)s 


where & is the volumetric strain rate caused by swelling and fi, fo, ... are predefined fields such as irradiation fluxes in 
cases involving nuclear radiation effects. Up to six predefined fields can be specified. 


Volumetric swelling cannot be used to define a rate-dependent gasket behavior. 


Input File Usage: 
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* SWELLING 
Abaqus/CAE Usage: 
Property module: material editor: Mechanical > Plasticity > Swelling 


Anisotropic swelling 


Anisotropy can easily be included in the swelling behavior. If anisotropic swelling behavior is defined, the anisotropic swelling 
strain rate is expressed as 


— 


- SW 


SW 2 SW 22 
é r) 


WwW 2 SW 
Eq = Ey +&q + &33 = (rr + 122 +133) 3 


: SW 
where € _ is the volumetric swelling strain rate that you define either directly (discussed above) or in user subroutine CREEP. 
The ratios 711, 722, and 733 are also user-defined. The directions of the components of the swelling strain rate are defined by 
the local material directions, which can be either user-defined or the default directions (see Orientations). 


Input File Usage: 
Use both of the following options: 


* SWELLING 
*RATIOS 


Abaqus/CAE Usage: 
Property module: material editor: Mechanical > Plasticity > Swelling: Suboptions > Ratios 


rt 





25 








User subroutine CREEP 


User subroutine CREEP provides a very general capability for implementing viscoplastic models such as creep and swelling 
models in which the strain rate potential can be written as a function of equivalent pressure stress, p; the Mises or Hill's 
equivalent deviatoric stress, q; and any number of solution-dependent state variables. Solution-dependent state variables are 
used in conjunction with the constitutive definition; their values evolve with the solution and can be defined in this subroutine. 
Examples are hardening variables associated with the model. 





The user subroutine can also be used to define very general rate- and time-dependent thickness-direction gasket behavior. 
When an even more general form is required for the strain rate potential, user subroutine umaT (User-defined mechanical 


material behavior) can be used. 
Input File Usage: 
Use one or both of the following options. Only the first option can be used to define gasket behavior. 


*CREEP, LAW=USER 
*SWELLING, LAW=USER 


Abaqus/CAE Usage: 
Use one or both of the following models. Only the first model can be used to define gasket behavior. 


Property module: material editor: 
Mechanical > Plasticity > Creep: Law: User defined 
Mechanical > Plasticity > Swelling: Law: User subroutine CREEP 


Removing creep effects in an analysis step 


You can specify that no creep (or viscoelastic) response can occur during certain analysis steps, even if creep (or viscoelastic) 
material properties have been defined. 


Input File Usage: 
Use one of the following options: 


*COUPLED TEMPERATURE-DISPLACEMENT, CREEP=NONE 
*SOILS, CONSOLIDATION, CREEP=NONE 


Abaqus/CAE Usage: 
Use one of the following options: 


Step module: Create Step: 
Coupled temp-displacement: toggle off Include creep/swelling/ 
viscoelastic behavior 


Soils: Pore fluid response: Transient consolidation: toggle off Include creep/swelling/viscoelastic behavior 


Fr 
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Integration 


Explicit integration, implicit integration, or both integration schemes can be used in a creep analysis, depending on the 
procedure used, the parameters specified for the procedure, the presence of plasticity, and whether or not geometric 
nonlinearity is requested. 


Application of explicit and implicit schemes 


Nonlinear creep problems are often solved efficiently by forward-difference integration of the inelastic strains (the “initial 
strain” method). This explicit method is computationally efficient because, unlike implicit methods, iteration is not required. 
Although this method is only conditionally stable, the numerical stability limit of the explicit operator is usually sufficiently 
large to allow the solution to be developed in a small number of time increments. 


Abaqus/Standard uses either an explicit or an implicit integration scheme or switches from explicit to implicit in the same step. 
These schemes are outlined first, followed by a description of which procedures use these integration schemes. 


1. Integration scheme 1: Starts with explicit integration and switches to implicit integration based on either stability or if 
plasticity is active. The stability limit used in explicit integration is discussed in the next section. 


2. Integration scheme 2: Starts with explicit integration and switches to implicit integration when plasticity is active. The 
stability criterion does not play a role here. 


3. Integration scheme 3: Always uses implicit integration. 


The use of the above integration schemes is determined by the procedure type, your choice of the integration type to be used, 
as well as whether or not geometric nonlinearity is requested. For quasi-static and coupled temperature-displacement 
procedures, if you do not choose an integration type, integration scheme 1 is used for a geometrically linear analysis and 
integration scheme 3 is used for a geometrically nonlinear analysis. You can force Abaqus/Standard to use explicit integration 
for creep and swelling effects in coupled temperature-displacement or quasi-static procedures, when plasticity is not active 
throughout the step (integration scheme 2). Explicit integration can be used regardless of whether or not geometric 
nonlinearity has been requested (see General and perturbation procedures). 


For a transient soils consolidation procedure, the implicit integration scheme (integration scheme 3) is always used, 
irrespective of whether a geometrically linear or nonlinear analysis is performed. 


Input File Usage: 
Use one of the following options to restrict Abaqus/Standard to using explicit integration: 


*VISCO, CREEP=EXPLICIT 
*COUPLED TEMPERATURE-DISPLACEMENT, CREEP=EXPLICIT 


Abaqus/CAE Usage: 
Use one of the following options to restrict Abaqus/Standard to using explicit integration: 


Step module: Create Step: 
Visco: Incrementation: Creep/swelling/ viscoelastic integration: 
Explicit 


Coupled temp-displacement: toggle on Include creep/swelling/ 
viscoelastic behavior: Incrementation: Creep/swelling/viscoelastic integration: Explicit 


Automatic monitoring of stability limit during explicit integration 
Abaqus/Standard monitors the stability limit automatically during explicit integration. If, at any point in the model, the creep 
strain increment (2,42) is larger than the total elastic strain, the problem will become unstable. Therefore, a stable time 


step, Ats, is calculated every increment by 


Ca 
At, = 0.5——, 


al 
t 


cr 
where € is the equivalent total elastic strain at time t, the beginning of the increment, and € 


ol is the equivalent creep 


t 
strain rate at time t. Furthermore, 


where Gi is the Mises stress at time t, and 
E=2(1+v) (n: D?:n) + 2.52, 
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where 
n=04q|,/00 
is the gradient of the deviatoric stress potential, 
D2 
is the elasticity matrix, and 
E 
is an effective elastic modulus—for isotropic elasticity E can be approximated by Young's modulus. 


At every increment for which explicit integration is performed, the stable time increment, At,, is compared to the critical time 
increment, At., which is calculated as follows: 


errtol 


Cr =Cr 


e leae © le 


At, = 


The quantity errto/ is an error tolerance that you define as discussed below. If At, is less than At,, At, is used as the time 


increment, which would mean that the stability criterion was limiting the size of the time step further than required by 
accuracy considerations. Abaqus/Standard will automatically switch to the backward difference operator (the implicit method, 
which is unconditionally stable) if At, is less than At, for nine consecutive increments, you have not restricted 


Abaqus/Standard to explicit integration as discussed above, and there is sufficient time left in the analysis (time left > 50At 
). The stiffness matrix will be reformed at every iteration if the implicit algorithm is used. 


Specifying the tolerance for automatic incrementation 


The integration tolerance must be chosen so that increments in stress, Ao, are calculated accurately. Consider a one- 
dimensional example. The stress increment, Ao, is 


Ag = EAc® = E (Ae — Ae”), 


where Ae, Ae, and Je" are the uniaxial elastic, total, and creep strain increments, respectively, and E is the elastic 
modulus. For Ag to be calculated accurately, the error in the creep strain increment, Aeéerr, must be small compared to Ae# 
; that is 

t | i 


Det ne er 


Measuring the error in Ae as 
A cr e*| = "| At 
Eerr Nee — ly 
leads to 


“cr 2cr el __ 
(é bares l. At « Ae* = —., or 


igs 
errto E " 


You define errto/ for the applicable procedure by choosing an acceptable stress error tolerance and dividing this by a typical 
elastic modulus; therefore, it should be a small fraction of the ratio of the typical stress and the effective elastic modulus in a 
problem. It is important to recognize that this approach for selecting a value for errto/ is often very conservative, and 
acceptable solutions can usually be obtained with higher values. 


Input File Usage: 
Use one of the following options: 


*VISCO, CETOL=errtol 
*COUPLED TEMPERATURE-DISPLACEMENT, CETOL=errtol 
*SOILS, CONSOLIDATION, CETOL=errtol 


Abaqus/CAE Usage: 


Use one of the following options: 
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Step module: Create Step: 
Visco: Incrementation: toggle on Creep/swelling/viscoelastic strain error 
tolerance, and enter a value 


Coupled temp-displacement: toggle on Include creep/swelling/ 
viscoelastic behavior: Incrementation: toggle on Creep/swelling/ 
viscoelastic strain error tolerance, and enter a value 


Soils: Pore fluid response: Transient consolidation: toggle on Include creep/swelling/viscoelastic behavior: 
Incrementation: toggle on Creep/swelling/viscoelastic strain error tolerance, and enter a value 


rt 


Loading control using creep strain rate 


In superplastic forming a controllable pressure is applied to deform a body. Superplastic materials can deform to very large 
strains, provided that the strain rates of the deformation are maintained within very tight tolerances. The objective of the 
superplastic analysis is to predict how the pressure must be controlled to form the component as fast as possible without 
exceeding a superplastic strain rate anywhere in the material. 


To achieve this using Abaqus/Standard, the controlling algorithm is as follows. During an increment Abaqus/Standard 
calculates 7max, the maximum value of the ratio of the equivalent creep strain rate to the target creep strain rate for any 


integration point in a specified element set. If Tmax is less than 0.2 or greater than 3.0 in a given increment, the increment is 
abandoned and restarted with the following load modifications: 


Tmax < 0.2 p = 2.0 pola, OF 
Tmax > 3.0 p= 0.5 pola; 


where p is the new load magnitude and p,jq is the old load magnitude. If 0.2 < r,,,, < 3.0, the increment is accepted; and 
at the beginning of the following time increment, the load magnitudes are modified as follows: 


0.2 <tmax <0.5 p=1.5p 54; 
0.5< tmx <0.8 p=1.2p 4; 
0.8 < tmx <15 p=pog; or 
1.5 < fie 33.0 p= pyq/1.2. 


When you activate the above algorithm, the loading in a creep and/or swelling problem can be controlled on the basis of the 
maximum equivalent creep strain rate found in a defined element set. As a minimum requirement, this method is used to 
define a target equivalent creep strain rate; however, if required, it can also be used to define the target creep strain rate as a 
function of equivalent creep strain (measured as log strain), temperature, and other predefined field variables. The creep 
strain dependency curve at each temperature must always start at zero equivalent creep strain. 


A solution-dependent amplitude is used to define the minimum and maximum limits of the loading (see Defining a solution- 
dependent amplitude for superplastic forming analysis). Any number or combination of loads can be used. The current value 
of Tmax iS available for output as discussed below. 


Input File Usage: 
Use all of the following options: 


*AMPLITUDE, NAME=name, DEFINITION=SOLUTION DEPENDENT 
*CLOAD, *DLOAD, *DSLOAD, and/or *BOUNDARY with 
AMPLITUDE=name 

*CREEP STRAIN RATE CONTROL, AMPLITUDE=name, ELSET=elset 


The *AMPLITUDE option must appear in the model definition portion of an input file, while the loading options 
(*CLOAD, *DLOAD, *DSLOAD, and *BOUNDARY) and the *CREEP STRAIN RATE CONTROL option should appear in 
each relevant step definition. 


Abaqus/CAE Usage: 
Creep strain rate control is not supported in Abaqus/CAE. 


Elements 


Rate-dependent plasticity (creep and swelling behavior) can be used with any continuum, shell, membrane, gasket, and beam 
element in Abaqus/Standard that has displacement degrees of freedom. Creep (but not swelling) can also be defined in the 
thickness direction of any gasket element in conjunction with the gasket behavior definition. 
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Output 


In addition to the standard output identifiers available in Abaqus/Standard (Abaqus/Standard output variable identifiers), the 
following variables relate directly to creep and swelling models: 


CEEQ 
t : : 
Equivalent creep strain, ho 4/ =e" > é"" dt. 
CESW 


Magnitude of swelling strain. 


The following output, which is relevant only for an analysis with creep strain rate loading control as discussed above, is printed 
at the beginning of an increment and is written automatically to the results file and output database file when any output to 
these files is requested: 


RATIO 

Maximum value of the ratio of the equivalent creep strain rate to the target creep strain rate, Tmax. 
AMPCU 

Current value of the solution-dependent amplitude. 
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